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ABSTRACT 



This thesis presents a method to calculate the power level 
of harmonics generated in DRFM due to sampling and 
quantization process. The Fourier series analysis is used for 
the calculation of the harmonic levels. The quantization 
process will generate the harmonics of the fundamental signal, 
which will be folded into the DRFM instantaneous bandwidth by 
the sampling process. Power level of the harmonics is 
dependent on the number of quantization bits of the analog- to- 
digital converter. Further, it is also dependent on the ratio 
of the signal frequency to sampling frequency. MATLAB 
programs for computation of harmonic power levels and plots of 
the harmonic power levels of the multibit DRFM are included in 
the thesis. Some applications of DRFM like deception jammer, 
broadband frequency source and radar simulator are also 
discussed . 
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I 



INTRODUCTION 



A. BACKGROUND 

The early generation radars were based on non-coherent 
detection techniques and therefore were susceptible to noise 
jamming. However, modern radars utilize coherent waveforms 
and sophisticated signal processing techniques such as pulse 
compression, doppler processing and digital correlators etc. 
These techniques provide additional processing gain to radar 
to discriminate unwanted signals including noise jamming. 
The jamming waveform can achieve the maximum processing gain 
if the waveform has similar characteristics as that of victim 
radar waveform. Comparable processing gain can be achieved by 
receiving, storing and retransmitting the victim radar's 
waveform with some suitable modulation. 

In the recent past, waveform storage was implemented with 
frequency memory loops (FML) . FML is an analog device using 
delay lines. Also, FML is generally not able to maintain 
coherence with the received waveform over a long storage 
period. Further FML is limited in terms of variety of 
deception techniques that could be implemented. 

The digital radio frequency memory (DRFM) is used for 
waveform storage. In DRFM, the incoming RF signal is sampled 
and stored in random access memory (RAM) . The RF signal is 
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regenerated before retransmitting to the victim radar. Since 
DRFM employs the digital techniques for waveform storage, it 
can maintain comparatively high coherence over a long storage 
period. Thus DRFM eliminates the signal fidelity 
deterioration with delay time. Because a complete pulse can 
be stored, pulse compression and phase coded signals with 
intra-pulse modulations can be stored and replicated with 
DRFM. When digital techniques are used, the flexibility to 
generate a variety of deception techniques is increased. Also, 
reprogrammability to meet scenario changes is an additional 
advantage . 

Some first generation DRFMs have been produced by 
several manufacturers. These first generation DRFMs are 
limited to single bit quantization. The single bit DRFMs have 
comparatively high harmonic levels. Thus useful jamming power 
is wasted in harmonics, reducing effective jamming power as 
well as providing a clue to the radar that it is being jammed. 
However, there are multi-bit, high performance DRFMs under 
development or in production. 

B. OBJECTIVE 

One disadvantage of DRFMs is that a sampling and 
quantization process generates the unwanted spurious/harmonics 
in the output. The objective of this thesis is to present a 
method for calculating the multi-bit DRFM harmonics/spurious 
signals. The operation of DRFM and its important 
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specifications are also discussed. Finally state of current 
technology, and some applications of the technology are given. 



C. RELATED WORK 

The information on DRFM technology is generally difficult 
to obtain. The reasons are 

• the technology is state-of-the-art and there is 
competition among manufacturers/contractors due to 
economics involved, 

• the technology is classified because of the military 
nature of the applications. 

Limited amount of literature is available in trade 
magazines and has been referred to in this thesis. 

D. OVERVIEW 

This thesis consists of five chapters and three 
appendices. Chapter I gives an introduction to the thesis. 
Chapter II describes the DRFM specifications. Chapter III 
provides the method for calculating the harmonics/spurious 
signals generated due to the sampling and quantization 
process. Chapter IV gives some of the applications for the 
DRFM technology, the current technology and specifications of 
a representative DRFM. Chapter V presents the conclusions and 
recommendations . 

Appendix A contains the MATLAB program for calculation of 
multi-bit quantized signal and plots of harmonic power levels. 
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Appendix B gives the MATLAB program for computation of 
harmonics of a squared error minimized quantized signal and 
the plots for the same. Finally Appendix C presents the 
MATLAB program for estimation of harmonic power levels of 
sampled and quantized signals and plots of harmonic power vs 
normalized input frequency. 
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II 



DRFM DESCRIPTION 



A. INTRODUCTION 

This chapter explains the functional organization of the 
DRFM and the major characteristics and requirements of the 
DRFM. 

B. DRFM OPERATION 

The simplified functional block diagram of the DRFM is 
shown in Figure 2.1. The RF channelizer downconverts the RF 
into a baseband signal for digitization. In single bit DRFMs 
the baseband signal is amplitude hard limited to obtain proper 
amplitude level for the digitizer. However, in multi-bit 
amplitude DRFMs, the baseband signal has to be processed 
through a complicated AGC circuit to control the proper 
amplitude levels required by the digitizer circuits. The 
multi-bit phase encoded DRFMs do not need the complicated AGC 
circuits because amplitude is not used for encoding. 

The sampler/digitizer is an analog-to-digital converter, 
which may encode either the amplitude or the phase. The clock 
frequency of the analog-to-digital converter is selected such 
that it is twice the maximum baseband frequency. The clock 
frequency is limited by the existing technology [Ref. 1] . 
Since memories cannot operate at the high clock frequencies, 
the multiplexer and demultiplexer circuits are required to 
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Figure 2.1 Functional organization of DRFM 



match the high data rate of the sampler to the read and write 
speed of the memory. The multiplexer circuits are 
parallel-to-serial converters and demultiplexer circuits are 
serial-to-parallel converters. The control circuits allow 
DRFM to generate the various deception waveforms. Control 
circuits also maintain the synchronous time reference for the 
proper operation of the DRFM as a whole. 

Since the basic function of the DRFM is to repeat the 
radar signals, the important requirement is to maintain 
sufficient signal fidelity at the output. Also DRFMs must be 
capable of meeting future needs for different deception 
techniques through software modifications. 

The following paragraphs discuss some of the major 
characteristics and requirements of DRFM. 

C. INSTANTANEOUS BANDWIDTH (IBW) 

The instantaneous bandwidth (IBW) is defined as the range 
of RF input signal frequencies at any one instant of time over 
which the DRFM system shall digitize, memorize, and recall the 
input signal without switching or tuning. The IBW of the DRFM 
determines the sampling clock rate of the system and also the 
memory size. 

Due to limitations of digital technology, the 
analog-to-digital converter and associated circuits can 
operate at maximum clock speeds on the order 1 GHz. Current 
DRFMs have an IBW on the order of 500 MHz. The IBW in the 
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range of 1 GHz - 2 GHz have been predicted for GaAs circuitry 
[Ref. 1]. 

D. FREQUENCY TUNING RANGE 

If the achievable IBW is less than the total bandwidth of 
expected threats, then the DRFM may have to be steered to the 
threat frequency. This requires the system to position the 
IBW anywhere in the desired threat spectrum using suitable 
local oscillators. Depending on the conversion hardware, it 
is feasible to cover an octave or possibly more than octave 
bandwidth. DRFM can be used to jam several threats 

simultaneously. Also, the unwanted signals in the IBW can be 
filtered out using a tunable bandstop filter. Signal gating 
in the time domain can also be used to select a threat from 
among several signals for storage in DRFM. 

E. LOCAL OSCILLATOR REQUIREMENTS 

The type of local oscillator required is dependent on 

■ rate at which IBW must be stepped and settled, 

■ the frequency accuracy of the DRFM output and 

■ maximum storage time for pulse 

The stepping of the IBW depends on the slew rate of the local 
oscillator and its settling time at the desired frequency. 
Similarly, frequency accuracy of the DRFM output is dependent 
on the short term frequency -stability. The maximum storage 
time of the pulse is dependent on the long term frequency 
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stability of the local oscillator. The type of local 
oscillators that meet the above requirements are fixed 
frequency oscillators like dielectric stabilized oscillator 
(DSO) , digitally tuned voltage controlled oscillator (VCO) , 
crystal oscillator-comb generator, or a frequency synthesizer. 

The highly accurate replicas are required for deception 
jamming of pulse doppler radar. This leads to the requirement 
that residual FM on the local oscillator should be low. Also, 
the FM on local oscillator should not generate the inter- 
modulation products in the required pulse doppler bandwidth of 
victim radar. A synthesized source appears to be good choice 
in this regard, but it comes at the expense of settling speed. 



F. DRFM SENSITIVITY 

As DRFM Utilizes down conversion hardware similar to the 
superheterodyne receiver, it is subjected to similar 
constraints involving sensitivity, bandwidth and false alarm. 
Proper RF and IF amplification are to be used to meet the 
required minimum signal-to-noise ratio. The sensitivity of 
the DRFM can be calculated by computing the signal-to-noise 



ratio (SNR) 


required 


to hold the 
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rates 


to 


reasonable 


values . 


Figure 2.2 


shows the 


plot 


of 



instantaneous bandwidth vs sensitivity for an SNR of 13 dB. 

As shown in the plot, the sensitivity is dependent not 
only on IBW, but also the noise figure of the system. So, it 
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In atari taneous BanciwIdtH vs DRF"M Sensitivity 




In sta nta nsous BandwIdtH 



Figure 2,2 DRFM sensitivity vs instantaneous bandwidth 

is important that RF-IF implementation incorporates low noise 
components to achieve the required sensitivity. 

G. DYNAMIC RANGE 

The dynamic range of the DRFM is dependent on the number 
of bits used in the quantization. Single bit DRFM can operate 
over a large dynamic range compared to multi-bit DRFM without 
additional circuitry. However, the multi-bit DRFMs can be 
made to operate over the required dynamic range with 
additional hardware, like instantaneous automatic gain control 
circuits, to provide proper input for the analog-to-digital 
converters. The input to the analog-to-digital converter has 
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to be maintained at the required level for spurious 
suppression at the output of the DRFM. 

H. SAMPLING REQUIREMENTS 

The sampling rate is the limiting factor in the design of 
DRFM. From the Nyquist sampling theorem, the signal has to 
be sampled at a frequency equal to or greater than twice the 
highest frequency component of the signal or twice the 
bandwidth of the signal. Hence in the case of DRFM, the 
sampling rate should be at least twice the IBW. If this 
requirement is met, it is ensured that proper reconstruction 
of the signal is possible. For example, if the IBW required 
is 500 MHz, the sampling rate should be greater than or equal 
to 1000 MHz. 

The IBW can be effectively increased by quadrature 
sampling [Ref. 2]. A block diagram of quadrature sampling is 
given in Figure 2.3. In quadrature sampling scheme, one 
sample per cycle is taken from both I and Q channels. This 
results in taking two samples per cycle using two analog-to- 
digital converters without increasing the sampling rate of 
analog-to-digital converters. Thus the sampling rate will be 
same as the IBW instead of twice the IBW. This is an 
important design consideration since the same sampling rate 
will achieve twice the IBW. 

However the IQ channel DRFM has some disadvantages. Due to 
the gain and phase mismatch between the I and Q channels the 
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Figure 2.3 Quadrature sampled DRFM block diagram 



output signal distortion may be increased. Another drawback 
is increased circuit complexity due to doubling of channels. 

I. QUANTIZATION REQUIREMENTS 

The sampler/ quantizer can be either an amplitude or a 
phase quantization analog-to-digital converter. The analog- 
to-digital converter should operate at a clock speed 
compatible with instantaneous bandwidth of the DRFM. Also the 
clock speed should be compatible with memory devices. 
Quantization is the most crucial function in terms of the 
achieving the waveform fidelity of the DRFM. The level of 
harmonic suppression required dictates the selection of the 
sampling rate and also the number of bits of quantization. 
For a given sampling rate, the harmonic signal level is 
reduced as the number of quantization bits are increased. 
Harmonic levels of multi-bit amplitude encoded DRFM are 
discussed in Chapter III. 

J . STORAGE REQUIREMENTS 

The memory size of the DRFM is determined by the 
instantaneous bandwidth and the maximum signal length to be 
stored (which is equal to longest pulse width) . The total 
number of bits to be stored (N) is given by 
N = f • N • N. 

s s b 

where 

fg = sampling rate = 2 • IBW 

Ng = maximum pulse width of the threat 

Njj = number of bits per sample 
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Assuming that one bit per sample is collected from the A/D and 
that the sampling rate is 50 MHz, one bit will be available at 
every 20 nanoseconds. If the longest pulse width is to be 
stored is 100 microseconds, the storage required will be 5000 
bits. 

Due to technology limitations, the IBW required to meet 
some types of threats is still difficult to achieve. However, 
multiplexing schemes can be used to provide higher IBW with 
achievable sampling rates. There are two multiplexing schemes 
available known as frequency multiplexing and time 
multiplexing. 

1. Frequency Multiplexing 

In frequency multiplexing, the signal is divided into 
a number of channels of reduced bandwidth. Each channel is 
sampled and converted to digital bit stream and stored. To 
reconstruct the signal, each channel is then recalled and an 
analog signal is reconstructed. By proper mixing and 
filtering, the composite threat signal is regenerated. 

2. Time multiplexing 

In time multiplexing, the signal is sampled with 
several time-phased clocks to produce several data streams 
which are then stored. The data streams are interleaved in 
an appropriate time sequence to reproduce the original 
signal. Duplex and quadraplex time multiplexing is possible 
to achieve higher sampling rates. The frequency multiplexing 
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is comparatively difficult process and hence the time 
multiplexing scheme is preferred. 

In this chapter, we have discussed the major 
characteristics and requirements of the DRFM. In the next 
chapter, we will calculate the levels of harmonic signals 
generated due to the quantization process. 
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SPURIOUS SIGNAL ANALYSIS 



A. INTRODUCTION 

The spurious signal is generated in the DRFM due to 
quantization and sampling. The spurious levels are of 
importance since they may create the recognizable signature of 
a jammer using DRFM. Also the power in spurious is not useful 
in the jamming. The theoretical analysis of the spurious 
signals generated in DRFM is elaborated in this chapter. 
Also a method to calculate the spurious signal level generated 
due to quantization is given. 

B. WAVEFORM SAMPLING AND QUANTIZATION 

The radar waveforms are continuous in amplitude and time. 
Therefore, the waveforms are to be converted to digital format 
for storing the radar waveform in the memory. The conversion 
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process involves two operations: sampling and quantization. 
In principle, operation of sampling and quantization can be 
performed in any order. However, in practice, conversion 
process is performed by the analog-to-digital converter and 
sampling operation is performed first. The analog-to-digital 
converter can be modeled as quantizer and sampler as shown in 
Figure 3.1. 

1. Quantization of the signal 

The quantization process can be considered to result 
from the application of signal to a staircase transducer. A 
staircase transducer converts the input voltage into binary 
levels representing the amplitude of the input signal. 
Assuming the input signal is sinusoid, the quantizer maps the 
amplitude of the sinusoidal signal to discrete levels 
represented by the binary number, i.e., the real valued 
instantaneous amplitude 'a' of the signal is mapped to a 
number k by the transducer. Number k is a binary number where 
0<k<N, N = 2'" where m is the number of bits of the 
analog-to-digital converter. Therefore the output of the 
quantizer is a discrete level signal instead of continuously 
varying signal. Because of the nonlinear characteristics of 
the quantization process, the harmonic distortion is expected 
in the output signal of the quantizer. Assuming that the 
sinusoid is quantized to one bit, that is to two levels 
representing positive and negative half cycles of the 
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sinusoid, the sinusoid input signal is converted to a square 
wave at the output of the quantizer. The harmonic content of 
the square wave signal is shown in Figure 3.2. It can be seen 
that the power of the each harmonic reduces monotonically . 

A method for calculating the harmonic content of the quantized 
signal is presented in the next section. 




Figure 3.2 Harmonics of the square wave 



2. Sampling of the signal 

The Nyquist sampling theorem states that an analog 
signal should be sampled at a rate greater than twice its 
highest frequency component. If this condition is satisfied, 
the continuous signal can be theoretically reconstructed 
completely from its samples as per sampling theorem. 

The sampling theorem, in its simple form, assumes 
infinite quantization. However, due to circuit limitations, 
the infinite quantization is physically not possible. So, the 
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harmonic distortion is inevitable in the DRFM sampling and 
quantization process. 



Quantized signal components 
Folded components due to fs 
Folded components due to 2fs 




Figure 3.3 Harmonics due to quantization and ideal sampling 



The sampling theorem also assumes that sampling width 
is infinitesimally small as represented by Dirac delta 
function 5(t) . As the sampling function is periodic in time, 
the input spectrum is replicated around the integer multiples 
of sampling frequency in the frequency domain. In Figure 3.3, 
the input signal spectrum has frequency components f^, ... 
generated by the quantization process. (Only up to 15th 
harmonic is shown in the figure for convenience) . The input 
spectrum is replicated around sampling frequency f^ as (fg'fo^ ' 
(f^-3f„) and so on, and also (f^+f^)/ (fc+3f.) and so on. 
similarly, the replication around twice the sampling frequency 
2fg will be (2fg-fjj) , (2fg-3f^) , so on. As the band of interest 
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is half the sampling frequency fg/2, it can be seen from 
Figure 3.3 that the harmonic signals (f^-lSf^j), (fg-13f^), 
(fg-llf^j), (fg-9fjj) are folded into the signal bandwidth due to 
sampling. Also, it is to be noted that the spectral lines due 
to integer multiples of sampling frequency 3f^, 4f^, . • • etc 
will also be folded into the signal bandwidth, but will be of 
much smaller magnitude. 

3. Effects due to finite width sampling 

In practice, the sampler is clocked with a periodic 
clock signal of finite width and hence the samples are of 
finite width repeating at the clock period of T^. The 
effect of quantization and finite width sampling on input 
sinusoid, quantized to single bit, is shown in Figure 3.4. As 
can be seen from the Figure 3.4, the harmonics of the 
fundamental generated by the quantization process are folded 
into the passband of DRFM (which is d.c. to half the sampling 
frequency) . 

As in case of infinitesimally small width sampling, 
the spectral lines generated by the integer multiples of the 
sampling frequency are also folded into the DRFM passband. 
However, the effect of finite width sampling will be to modify 
the spectral levels by the sine function as shown in the 
Figure 3.4. The harmonic distortion will be proportional to 
ratio of the signal frequency to sampling frequency. If the 
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ratio of the signal frequency to sampling frequency is low, 
harmonic distortion will also be low. 

Though the input signal is band-limited to frequency less 
than half the sampling frequency, the quantization process 
generates the harmonics of the fundamental signal and the 
sampling process folds them into the passband. If the input 
signal has wider spectrum, it is to be expected that the 
quantization process will generate harmonics of each of the 
spectral line in the input signal. Further, the sampling 
process will fold the intermodulation products of harmonics 
and sampling frequency into the instantaneous bandwidth of the 
DRFM. As the sampling and quantization processes are 
essential, the requirement will be to reduce the harmonic 
levels due to quantization such that overall spurious outputs 
of the DRFM is with in acceptable levels. 

C. COMPUTATION OF HARMONICS DUE TO QUANTIZER 

This section presents a method to calculate the harmonic 
signal levels generated by the quantizer. 

In order to calculate the harmonic distortion introduced 
by the quantizing process, we consider a sinusoidal signal 
x(t) = AjjSinOjjt . A quantized sinusoidal signal can be 
represented by the quantized signal Xq(t) as shown 
Figure 3.5. If m is number of bits of quantization , then 
each half cycle of the sinusoid is quantized into N steps 
where N = 2’"'^. The quantized signal x^(t) can be considered 
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as superposition of N pairs of rectangular step signals of 
period equal to the period of the sinusoid. The amplitude 
of each rectangular step is A^N. The width of step k is t^, 
in each half cycle, and is given by 

sin~^ ...(3.1) 

^ 2 N 




Figure 3.5 Quantized signal step pairs 



The equation 3.1 is equivalent to 



tu = 



— sin'^ — 



0 ), 



N 






2 

2nf^ N 

lo sin-^A]SzlL 

TZ N 
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For two bit amplitude quantization, number of pairs N 
amplitude step = A^2 and k = 1,2. Then for N = 2, k =1 
width t|^ is given by 



tjc 




The waveform for N = 2, k = 1 is as in Figure 3.6 









and for N = 2, k = 2, the width t,^ is given by 




3 



= 2 , 
, the 
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The waveform for N = 2, k = 2 is as in Figure 3.7 



T T 

K 3 ^12 




The quantized waveform of the 2 bit amplitude quantized 
sinusoid is given in Figure 3.8. 



To 




Figure 3.8 Quantized combined waveform for N = 2, k = 1,2 
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The k'th rectangular step pair may be represented by the 
Fourier series as 



where is given by 

(,) — -ACjt 

■®njc =^ — -I sinnw^tdC n = l,3,5,7 

N n Jacu ° 



where At,^, offset of the of k'th step is 



A tj, = J^sin-^ ( JJLJA. ) 
* co„ N 



(3.2) 



...(3.3) 



. . . (3.4) 



Then B„ is 



B 



4A 



uN{2n+l) 



cos [ (2n+l) sin'^ ( ) ] 

N 



(3.5) 



Hence the Fourier series of the k'th step pair is 

=E“o^^77l^-rrC°s[(2n + l) sin-M-^^^)] 

z^n=o 7iiy(2n + l) N 

• sin (2n+l) o)(jt ...(3.6) 

For N pair of rectangular pulses, we obtain Fourier series of 
the quantized signal Xq(t) as 

E n 4^ 



-cos[(2n+l) sin~^( )] 

Ki\r(2n+l) N 



sin (2n+l) co„t 



= ELo^2n*iSin(2n+l)<o^t 



. . . (3.7) 
. . . (3.8) 
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where is the amplitude of the (2n+l) 'th harmonic and is 

given by 



^2n + l' 



— rYl cos [ (2n+l) sin'M-- --^^ ) ] 

m7\^(2n + l) ^^=1 N 



(3.9) 



The signal power in each of the harmonics, given by 
(^2n+i)^' can be calculated using equation 3.9. The harmonic 
power for each harmonic is shown in Appendix A in Figures A.l 
through A. 18. It can be seen from the figures that 

• the third harmonic is predominant 

• the total harmonic power decreases as the number of bits 
of quantization is increased. 

Table 3.1 shows the maximum harmonic power in dB relative to 
fundamental and total harmonic power relative to fundamental 
for different quantization levels. 

As the third harmonic is predominant and it will be folded 
back into the IBW due to sampling, it is advisable to move 
the strongest harmonic from the third to a higher harmonic. 

D. REDUCTION OF HARMONIC POWER 

As discussed in the last section, the harmonic power level 
in the third harmonic is maximum for all the quantization bit 
levels. Also this harmonic will be folded into the passband 
by the sampling process. Hence, it is necessary that the 
harmonic signal with maximum harmonic power level be as far 
from the fundamental as possible so as not to fold into the 
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TABLE 3.1 

COMPARISON OF QUANTIZATION BITS AND HARMONIC POWER LEVEL 



Number 
of bits 


Third harmonic 
power level 


Total harmonic power 
in 99 harmonics 


1 


- 9.5 dBr 


- 6.4 dBr 


2 


- 14.9 dBr 


- 12.1 dBr 


3 


- 21.3 dBr 


- 18.0 dBr 


4 


- 28.9 dBr 


- 24.1 dBr 


f 

5 


- 34.7 dBr 


- 30.6 dBr 


6 


- 41.3 dBr 


- 38.2 dBr 



passband. The following paragraphs present certain methods 
to reduce the harmonic power in the output of the DRFM. 

1. Harmonic reduction in single bit quantization DRFMs 
The spur reduction in single bit DRFM is based on the 
relationship between the level of odd harmonics of the 
rectangular wave and its duty cycle [Ref. 3], By adjusting 
the duty cycle to 1/N, the harmonic spur N can eliminated at 
the expense of generating even harmonics. If the duty cycle 
is varied from less than 0.5 to greater than 0.5, the energy 
is transferred from odd harmonics to even harmonics, but due 
to reverse duty cycle, the even harmonics are phased out. 

It is also possible to reduce the harmonics by use of 
controlled time jitter in the sampling clock [Ref. 4]. The 
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time jitter t of the clock should be controlled such that t„ 
is an identically distributed independent random variable with 
a uniform probability density function over [-nT^2, nT^2], 
where 0 < n <= 1, and is the time period of sampled signal. 
It is stated that power level of the fundamental signal is 
also reduced to some extent. It is to be noted that this 
method can also be used in case of multi-bit quantization. 

2. Harmonic reduction multi-bit quantization DRFMS 

Consider the input sinusoid signal quantized to 2'"'^ 
amplitude levels where m is the number bits of the analog to 
digital converter in each of the half cycles as shown in 
Figure 3.9. The amplitude step of quantization is 'a' and 
the transition angle at the quantization point is 0 as shown. 




Figure 3.9 The quantization levels and angles of signal 
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Assuming positive and negative half cycle symmetry in 
the quantized signal, the Fourier series of the 2’"'^ level 
quantized signal is given by 

=5^“^^i3,,sinno)„t ...(3.10) 

where is given by 

2 r®2 

b„ = — [/ asinnxdx+l 2asinnxdx+ 3asinnxdx+... 

" TC Jo J02 J03 

+ ^ 2'""^a sin nxdx + ... + f a sinnxdx] . . . (3.11) 

*'71—02 

i . e. , 

b^= .. . — ( 1 -cos ( ( 2 n+l) k) ) 

^ ( 2 n+l) 7 t 

• [1 + cos ( 2 n+l) 02 + cos ( 2 n+l) 63 

+ . . . cos (2n+l) 02»-i] . . . (3.12) 



For two bit DRFM, the harmonic amplitude levels are given by 

[1 -^cos(2n+l)02] ...(3.13) 

" ( 2 n+l)Ti ^ 

For three bit DRFM, the harmonic amplitude levels are given by 
b„~— — ^-^-r — [1 + cos (2n+l) 0, + cos ( 2 n+l) 0,+cos ( 2 n+l) 0.] ...(3.14) 

(2n+l)7t 2 3 4 

The step amplitude level 'a' and transition angle 0^, 
s = 1 , 2 , .. . can be selected such that staircase output is best 
fit to the input sinusoid. To determine the best fit step 
amplitude 'a' and transition angles 0 s, we proceed as follows. 
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As shown in Figure 3.9, 
is as given below. 



the error in the different intervals 



Transition angle 
interval 



Error in the 
interval 



0 < 0 < 0 ^ 

0 . < 0 < 0 , 

03 < 0 < 0 ^ 



a - sin 0 
2 a - sin 0 
3a - sin 0 



Bj.., < e < ir/2 



2 "'-la - sin 0 



So the squared error over the complete cycle is given by 



input signal, the partial derivatives of the above equation 
with respect to 02 , 03 ,... 0 n,_i, and step amplitude 'a' must be 
equal to zero. Solving for step amplitude 'a' and 0j, 03 ,... 
0 ^.^ that meets the above requirement, we get following 
relations 



^ e^=4 f ^ (a-sin0) ^ d0 + (2a-sin0) ^ d0 + . . . 




7C 




For minimum error in the output signal with respect to 



sin 02 = 3a/2 
sin 03 = 5a/2 



sin 02 „,.i = ( 2 -"-l)a /2 



(3 . 16) 



and the step amplitude 'a' as below. 



Number of Bits 



Amplitude step 



1 

2 

3 

4 

5 



0.6366198 

0.4188161 

0.2332685 

0.1214300 

0.0616995 
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The harmonic signal levels can be calculated based on 
the above values for angle 0 s hnd amplitude 'a' in the 
equation 3.12. For two bit DRFM, the equation 3.12 reduces to 

[ 1 +cos ( 1262 ) ] 

” nit 

= [l+cos(38.92*n) ] ...(3.17) 

n 

The harmonic power in the higher bit DRFM can be 
calculated accordingly by using the above values in equation 
3.12. The plot of harmonic power in multi bit DRFMs 
calculated as per equation 3.12 is shown in Figures B.l 
through B .8 in Appendix B. However, it is to be noted that the 
harmonic power as shown is dependent on the sampling point in 
the RF cycle and amplitude at the sampling point. In a high 
frequency sampling system, the point of sampling becomes 

critical and hence calls for a high stability clock for 

sampling. The plots in Figures B.l through B .8 show that 

■ the third harmonic is not the worst offender 

■ the harmonic power is not decreasing with 

increasing harmonic number 

Table 3.2 gives the harmonic power level and total 
harmonic power for versus number of bits. 
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The results in the Table 3.2 are more conservative because the 



following factors contribute towards reduction of the harmonic 
content in actual design [Ref. 5] . 



• rounding of the quantized waveform due to stray 
capacitances in the circuit. 

» generation of even harmonics due to nonsymmetric quantized 
waveform will reduce the odd harmonic power, 

♦ drift in the L.O. frequency spreads the harmonics over the 
IBW. 



TABLE 3.2 

NUMBER OF QUANTIZATION BITS AND HARMONIC POWER 



Number 


Worst 


Power level of 


Total Harmonic 


of bits 


Harmonic 


worst harmonic 


power w.r.t 
fundamental 


i 2 


9 


- 18.1 dBr 


5.44 % 


3 


19 


- 25.5 dBr 


1.17 % 


4 


43 


- 34.3 dBr 


0.27 % 


5 


45 


- 45.9 dBr 


0.63 % 



Also, note that the results in the Table 3.2 are more 
optimistic than actual as the harmonic power level dependency 
on the relation between signal frequency and sampling 
frequency is not considered. 

In the next section, we discuss the dependence of 
harmonic power level with respect to signal frequency and 
sampling frequency. 
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E. EXPRESSION FOR SAMPLED AND QUANTIZED SIGNAL 

The quantized signal x^(t) is sampled with the 
clock signal x^(t). Thus the sampled quantized signal Xg(t) 
is given by 

Xs(t) = Xq(t) • xjt) ...(3.18) 

where Xq(t) is the quantized signal as derived before. The 
clock signal x^(t) is a signal as shown in Figure 3.10. 












Ac 




-Te 


0 




Tc 



Figure 3.10 Clock signal 



The Fourier series of the clock signal x (t) is given by 



i7?C0j,X^ 



[1-h2V“ ^ 






• cosnco„t] 



....(3.19) 



Replacing clock signal x^.(t) and quantized signal x^(t) with 
their Fourier series expansions, we get x^(t) as 



X (t) =A [1+2V“ 

S C rp ^m -1 

c 



Sin 



;t70)^Xc 






4A, 



-cos/no>„t] 



cos [ (2/3+1) sin'M-^^y^) ] 
■KN{2n+l) N 



sin (2/3+1) o>„t] 



(3.20) 
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By expansion and trigonometric manipulation we have 



E! — ff£_^cos[(2n+l) sin-1 

s Z-^/ 3 t=o - X 



TiN(2n + l) 

• sin (2n+l) 



N 



E » 

m=l 2L/n^o 2-^k=l 



Sin 






T. TtN{2n+l) /no)_T, 



cos [ (2n+l) sin~i ( )3 sin[/no)^t+ (2n+l)w^t] 



E og C 

^n=0 



Sin 



inoo^T^ 



K//(2n+l) 

/ u_'| \ 2 

• cos [ (2n+l) sin ^ ^ ^ ^ ) ] sin [/nw^t- (2n+l) co^t] 

... ( 3 . 21 ) 



The first term in the above equation 



E « 

n=0 2^k=l 



4^. 



r. 7rT/(2n+l) 



cos 



[ (2n+l) sin'i { ) ] 



N 



* sin [ (2n+l) cOq t] ...(3.22) 

gives all the harmonics of the fundamental frequency o^. The 

harmonic frequencies beyond the filter bandwidth f<-/2 are 

rejected by the filtering process. Hence only the frequency 

terms (2n+l)(i)^ less than 0^/2 are to be considered. 

The second term in the above equation 



E oo C 

in=l 2^n=0 2^k-l 



4A, 



Sin 






T„ 'KN{2n+l) 



-COS 



[(2n+l)sin-i(-%^)] 



• sin [irjco^t+ (2n+l) co^t] 



N 

. . . ( 3 . 23 ) 
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contribute 



does 



not 



as higher order frequencies 



[(2n+l)o) + mo) ] generated by this term are rejected by the 



filtering process because they are above o)(,/2 
term in equation (3.21) can be neglected. 

The third term in equation (3.21) is 



So the second 



^=1 2^k-l 



4A. 



Sin 






nN(2n+l) 






■cos [ ( 2 I 2 ■^ 1 ) sin"^ ( ) ] 

N 



sin[;nG)^t- (2n+l) u>^t] 



. . (3.24) 



and it generates the inter-modulation products as in a mixing 
operation and contributes to the harmonic power if the 
[2(n+l)o)^ - mo)^] is less than o)^/2. 

Finally, the sampled signal x^(t) can be written as 



xAt) =y‘ 

^ Z^n=0 



r_ nN{2n+l) 



cos 



[ (2/2+1) sin-^-^^1^] 



N 



•sin ( 2 / 2 + 1 ) G)-C 



E <» ^^00 

in=l 2^k^\ 






4A^ 



Sin 



/7?C0^T^ 



T- nN(2n+l) /no)_t^ 



■cos [ ( 2 / 2 + 1 ) sin ^ ( ) ] • sin [/nco^C - ( 2 / 2 + 1 ) co^t] 

...(3.25) 



The equation (3.25) shows that the inter-modulation 
products of the harmonics of the signal frequency and 
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harmonics of the clock frequency will contribute to the spurs. 
It is convenient to rewrite the equation (3.25) as 



tZ-O tH-l - (2n + l) G)„t] 



. . (3.26) 



where ^ is defined as 






4A, 



■ n,k 



, -cos [ (2n+l) sin'^ ] 

nN{2n+l) N 



...(3.27) 



and Q i, is defined as 

n, k 



^m,n, k 






4 ^. 



Sin 






TiN[2nkl) 






- — cos [ (2n+l) sin-^-^^-ii] ...(3.28) 

N 



The plots of harmonic power in dB vs frequency as per 
equation 3.26 is given in Appendix C. The Figures C.l through 
C.16 show the plot for 1, 2, 3 and 4 bit amplitude quantized 
sampled signals. As can be seen from the plots, the harmonic 
power folded into the instantaneous bandwidth is dependent on 
the relation between input signal and the sampling frequency. 
If the ratio of the signal frequency to sampling frequency is 
low, the harmonic power relative to the fundamental is also 
low. Table 3.3 shows the relationship of normalized input 
frequency to harmonic power for 1 through 4 bits sampled 
quantized signal. 

As can be seen from Table 3.3, the third harmonic level of 
quantized signal will be approximately equal to quantized 
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sampled signal only if the sampling frequency is high compared 
to signal frequency. As the ratio of signal frequency to 
sampling frequency increases, the harmonic power also 
increases. Tables C.l through C.4 list the harmonic power of 
1 through 4 bit sampled quantized signals. 

TABLE 3.3 

HARMONIC POWER LEVEL OF QUANTIZED AND SAMPLED SIGNAL 



1 Normalized 
Input 
Frequency 


Third Harmonic power level (dB) , 
Quantization Bits 


1 


2 


3 


4 ' 


Quantized 

only 


- 9.54 


- 14.96 


- 21.33 


- 28.00 


0.005 


- 9.45 


- 14.76 


- 20.86 


- 26.98 


1 

0.025 


- 8.86 


- 13.44 


- 18.28 


- 22.66 


0.05 


- 8.23 


- 12.18 


- 16.34 


- 19.10 


0.1 


- 7.01 


- 9.94 


- 12.96 


- 15.23 
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F. EXPRESSION FOR DRFM OUTPUT 



The stored signal is repeated with pulse repetition period 
Tj, (PRF of f^) , pulse width of 7 ^ and amplitude of A^. The 
Fourier series of such a waveform is 



X 



At) 



sin 



Jo) T 

r r 



^ 2 



■cosico^t ] 



...(3.29) 



So, the output of the DRFM is 

x^(t) = x^Ct) • x^(t) . . . (3.30) 

Replacing the Fourier series expansion for Xg(t) and x^(t) we 
have X|j(t) as 



^:,,jcSin(2n+l)a)„t 

" (2n + l) 0 )„t] ] 



T’ 2 ^ 1=1 r 



sin 



io) T 

■J-W w ^ 






■cos ( ] ...(3.31) 



2 
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By expansion and trigonometric manipulations 



(2n+l)G)^t) 



^EliELoK=i^->".Jc^sin(;7)0),- (2nH-l)0)J t 



Sin 



itO^T^ 



+Y^“ Y^“ p 2 

2-ri=i 2^n=0 2^k^l rp l(x) X 



• sin ( (2n-i-l) o>^+io)^) t 



sin 



io> T 

r r 



yv/ p 2 

Z^i*l Z^n^Q Z^k^i ^n,k rp Jg) T 



sin ( (2n+l) co^-lo^) t 



sin 



Jo> T 

-L % j. 



E oo y^oo y>(» y^vT/ _ z 2 

1-1 2^m-l 2^n-0 2^k-l ^^>ri,k rp lid ^ 



•sin (/nco^ - (2n+l) o>^ - t 



E o© y^oo ^ '^z'^Z 

1=1 2^m=l 2^n=0 2^k=i rp 



71 ^ sin i- 



io.'C, 



•sin (/ntOp - (2n+l) o)^ + J<o^) t 



. . . (3.33) 



The above equation defines the output of the DRFM completely. 
As can be seen from the equation, the effect of quantization 
and sampling the band-limited signal is to fold the harmonics 
produced by the quantization process into the signal passband. 
The harmonics outside the instantaneous bandwidth can be 
removed by placing an appropriate bandpass filter after the 
conversion to analog signal. But the harmonics folded into the 
instantaneous bandwidth of the DRFM can not be filtered out 
easily. 
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IV 



DRFM APPLICATIONS 



This chapter describes some of the applications of DRFM 
such as deception jammer, broadband frequency source, signal 
analyzer and simulator. Some ECM techniques using DRFM are 
also described here. 

A. ECM SYSTEM USING DRFM 

This section will describe a DRFM based system which will 
capture radar signals in a band covering 8 to 16 GHz. The 
block schematic of the system is shown in Figure 4.1. The 
DRFM used in this system has an instantaneous bandwidth of 500 
MHz. As the frequency range required by the system is 8 GHz, 
it is necessary to steer the instantaneous band (500 MHz) of 
DRFM to the threat signal frequency. Instantaneous frequency 
measuring (IFM) receiver is used for measuring the incoming 
signal frequency. The measured frequency is used to chose a 
local oscillator from the local oscillator bank. The local 
oscillator bank contains a set of fixed frequency oscillators 
as given in Table 4.1. The delay line in front of the DRFM 
will provide necessary time delay required for setting the 
local oscillator frequency. 

The incoming signal is first down-converted to a 500 MHz 
band centered at 1 GHz by selecting an appropriate local 
oscillator. The resulting down-converted signal is quadrature 
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RX ANTENNA 
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Figure 4.1 Frequency Steered DRFM block schematic 



mixed against 1 GHz local oscillator to provide +/- 250 MHz 



TABLE 4.1 SELECTION LIST OF LOCAL OSCILLATOR 
FREQUENCY, IF FREQUENCY BAND FOR EACH SUB-BAND 



Freq range 


L.O. 


IF 


(GHz) 


(GHz) 


(GHz) 


8. 0-8. 5 


9.25 


1.25 -0.75 


8. 5-9.0 


9.75 


1.25 -0.75 


9. 0-9. 5 


10.25 


1.25 -0.75 


9.5-10.0 


10.75 


1.25 -0.75 


10.0-10.5 


11.25 


1.25 -0.75 


1$. 5-14.0 


11.75 


1.25 -0.75 


11.0-11.5 


11.25 


0.75-1.25 


11.5-12 . 0 


11.75 


0.75-1.25 


12.0-12.5 


14.25 


0.75-1.25 


12.5-13.0 


13.75 


1.25 -0.75 


13.0-13.5 


14.25 


1.25 -0.75 


13.5-14.0 


10.75 


1.25 -0.75 


14.0-14.5 


10.25 


0.75-1.25 


14 . 5-15.0 


13.75 


0.75-1.25 


15.0-15.5 


14.25 


0.75-1.25 


15.5-16.0 


14.75 


0.75-1.25 



baseband in-phase (I) and quadrature (Q) signals. The I and 
Q signals are passed through signal conditioner circuits. 
Each of the conditioned I and Q signals are sampled and 
quantized with 500 MHz clock in analog-to-digital converter 
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(A/D) . The samples are then stored in a random access memory 
(RAM) for later signal reconstruction. To replicate the 
signal, the inverse process is implemented. The stored 
digital samples are clocked out at 500 MHz rate to the 
digital-to-analog converter (D/A) for generation of an analog 
baseband signal. The regenerated baseband signal is 
upconverted with the same local oscillators which were used in 
the down conversion process. Thus DRFM will replicate the 
instantaneous frequency of the received signal at the output. 

B. ECM TECHNIQUES 

Basic ECM techniques such as range gate stealing, velocity 
gate stealing, cover pulse noise jamming, false target 
generation, etc., can be implemented with DRFM. 

1. Range Gate Stealer (RGS) 

Range gate stealing is an effective ECM technique 
against the tracking radars. The range deception is 
accomplished by transmitting the stored signal to capture the 
radar automatic gain control circuits (AGC) . After the 
capture of the AGC, the stored signal is transmitted with a 
successively increased delay to generate the required range 
gate pull-off. At the end of maximum delay, the cycle of range 
gate stealing is repeated. The increment of range delay can 
be made to imitate any pattern required to simulate different 
types of platform maneuvers. 
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Due to storage of complete pulse width by the DRFM, 
the phase coded and chirped waveforms can be coherently 
reproduced with controlled pulse width and phase. Due to long 
storage times available, it is feasible to generate range gate 
pull-in for generating approaching targets. 

2. Velocity Gate stealer (VGS) 

The velocity deception or velocity gate stealing is 
effective ECM against the Doppler sensing radars like CW radar 
and pulse Doppler radars. The VGS is accomplished by 
repeating a frequency shifted replica of the victim radar's 
signal. The frequency shift is initially programmed so that 
repeated signal is with-in the passband of the Doppler filter 
containing the target return. This allows the jammer to 
capture the radar's AGC circuits. The repeated signal is 
then slowly walked off in frequency to the expected maximum 
Doppler frequency of the radar. The repeated signal is then 
removed, forcing the radar to reacquire the target. 

In DRFM, the frequency of the signal can be shifted, 
by introducing the required frequency change in the 
upconverting local oscillator with respect to down converting 
local oscillator. The Doppler frequency added to the 
upconverting local oscillator can be generated digitally to 
the required accuracy. As the stability of the of DRFM output 
is determined by the stability of the LO and sampling clock, 
short term stability of few Hz can be easily achieved. The 
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multi-bit DRFMs become a necessity if the system is to be used 
for velocity gate stealing since the harmonic signal and 
spurious signal levels decrease as the number of quantization 
bits are increased. 

3. Coordinated RGS/VGS 

RGS or VGS, by themselves, may not deceive the pulse 
Doppler radar. The pulse Doppler radar may check the 
consistency of target velocity obtained by differentiating 
range data with that obtained from target Doppler data. 
Hence, to deceive pulse Doppler radar, RGS and VGS are to be 
coordinated so as to provide a realistic target motion. 

The ability of DRFM to coordinate is excellent as 
digital techniques are used for signal storage and 
reconstruction. Precise RGS is achieved due to timed readout 
of memory. VGS is coordinated by selecting an appropriate 
Doppler corresponding to the RGS selected. Thus it will be 
difficult for pulse Doppler radar to reject the false target. 

4. Multiple False Target Generation 

As the readout from the memory is nondestructive, it 
is possible to regenerate any number of false targets. The 
false targets can also be generated at any range since the 
storage time is unlimited. The memory available can be 
segmented for use to generate different false targets. As the 
multiple targets are generated sequentially, it is possible to 
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toggle between different false targets in a coordinated 
fashion. However, care should be taken not to exceed the 
duty cycle limitations of the high power section. 

C. DRFM AS A FREQUENCY SOURCE 

The wide band signals required in radar systems can be 
generated by DRFMs. Digital data, required by DRFM, can be 
either from computer simulated data or from a real world 
digitized signal. Wide bandwidth signals such as phase coded 
or frequency chirped can be easily generated by this process. 
In radar designs, DRFM technology has been used to enhance the 
function of exciter of the radar by providing ready made and 
more versatile waveforms for various specialized functions 
[Ref. 6] . The main advantage of using DRFM technology will be 
that the signal waveform can be easily manipulated by changing 
the digital data in the DRFM memory. 

D. SIGNAL ANALYSIS USING DRFM 

The DRFM technology can also be used in receiving systems. 
The digital techniques are available for determining signal 
characteristics such as amplitude, phase and frequency of the 
signal. The block diagram of the digital receiver is given 
in Figure 4.2. 

Unlike conventional receivers, digital receivers can store 
the digitized signal information for long periods of time 
without degradation. The digital demodulation techniques can 
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Figure 4.2 Digital Receiver 

be applied to the stored signal to obtain information like 

amplitude, phase and frequency [Ref. 7]. If the incoming 

signal a(t)cos[w^t + p(t)] is quadrature sampled, it can be 

represented as 

I(t) = a(t) cos(p(t)) 

Q(t) = a(t) sin(p(t) ) 

The amplitude of the signal can be obtained by 

a ( t) 

The phase information of the signal can be obtained by 



<J)( t) = tan'M ] 



= tan'^ [ 



I{t) 

a{ t) sin(p( t) ) j 
a ( t) cos {pit)) ^ 



The p(t) contains both the carrier phase offset and the 
modulation frequency phase slope. The other signals, like FM 
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or phase coded, can be processed in a similar fashion as the 
pulsed signals. However, due to technology limitations, the 
real time analysis of the signals is not yet feasible over 
very wide bandwidths. 

Also the input data stored in the memory could be 
processed several times to obtain different signal 
information. For example, if a data set contains both strong 
and weak signals, it is possible to obtain information on a 
strong signal and then subtracting strong signal data from the 
stored data, the weak signal can be detected. This is useful 
in improving the dynamic range of the receiver. 

E. DRFM AS SIMULATOR 

The DRFM has been used in simulation for testing of radar. 
The DRFM can be used for simulation of targets, clutter, and 
ECM waveforms. The main advantages in this type of 
application are 

a) The digitized real world signal or the digital 
words of signal waveform simulated on computer 
can be stored in DRFM to generate the required 
signal . 

b) The system is generic in application; a single 
system can be used against different radars by 

modifying the data in DRFM memory. 

c) Due to digital techniques involved, the 
simulations are repeatable to high accuracy. 

d) Only a short reprogramming time is required to 
react to changes in scenarios. 
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Such systems have been integrated into certain ground based 
fire control radars for test purposes. 



F. SPECIFICATIONS OF DRFM 



Following is the typical 
DRFM [Ref. 8]. 

1 Frequency of operation 

2 . Instantaneous Bandwidth 

3. Dynamic Range 

4 . Memory length 

5. Output frequency 

6. Output pulse width 

7. Output power levels 

8. Operating Modes 
a . RGPO 

Dwell 

Walk off time 
Delay range 
Minimum delay 
Resolution 
Sweep 
Hold mode 



specification of a representative 

5.4 - 5.8 GHz 
400 MHz 
3 0 dB 

100 ns - 27 micro seconds 
Input +/- 1 KHz 
Input +/- 40 ns 
- 12 +/- 4 dBm 

0.5 - 5.0 seconds 
0.5 - 20.0 seconds 
1-30 micro seconds 
60 nanoseconds 
40 nanoseconds 
parabolic 
1-15 seconds 



b. False Targets (Preset replicas) 



Number of false targets 

Continuous wave (CW) 
External trigger 



2, 4, 8, 16 

spaced 26 microseconds 
apart 

continuous readout of 
stored pulse 
1 replica per trigger 
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G. TECHNOLOGY TRENDS 

The DRFM technology has developed in the past few years as 
a versatile means for effective ECM systems against coherent 
radars. DRFMs are being currently produced by several 
companies including Whittaker and Raytheon. Table 4.2 gives 
the list of DRFMs manufactured by different vendors with their 
brief specifications [Ref. 1]. 

The future developments in very high speed digital 
integrated circuits, analog-to-digital converters and 
digital-to-analog converters will further improve the 
performance of the DRFMs. The future trend appears to be the 
use of MICs and ASICs to reduce the weight and size the DRFMs. 
ITT Avionics has developed a DRFM that consumes 15 watts of 
power and is only 20 cubic inches in size. This is an 
improvement over the earlier 7 0 watt and 7 0 cubic inches 
version . 

With high speed intelligent processing built into the 
DRFMs, it will be possible to generate complex signals 
required for ECM applications. Efforts to design and 
fabricate GaAs/ECL DRFMs for higher circuit integration are 
reported. The sampling speed appears to be the significant 
factor driving the DRFM technology. DRFMs with sampling clock 
rates better than 1 GHz and IBW better than 1 GHz have been 
reported [Ref. 1]. 

It has been reported that the Air Force Institute of 
Technology (AFIT) and Wright Laboratory have plan to build a 
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DRFM onto a chip [Ref. 9] . The AFIT has used Hilbert transform 
to accomplish the amplitude and phase modulation. Both 
amplitude and phase of the signal is controlled on a sample by 
sample basis by this technique. The ultimate goal is to place 
all of the required converters, memory and control circuits on 
a chip. As noted earlier, the capabilities of the DRFMs will 
be enhanced further because of these efforts. 



TABLE 4.2 

LIST OF DRFMS WITH THEIR BRIEF SPECIFICATIONS 



Manufacturer 
and Model 
• Number 


RF Freq 
Range 
(GHz) and 
(IBW) 


A/D Bits 
and type 


Memory 

length 

in 

(micro 

sec) 


RF Output 


Power out 
dBm 


Modulation 

type 


Delay available 
(min - max) 


Anaren 

45010 


7.0 - 17.0 
(500) 


3, 

Phase 


200 


+4 


Doppler, 
Ampl i tude, 
Frequency 


15 ns - 3.07 ms 


CAL TCM-875 


7.0 - 11.0 
(1000) 


4, 

Ampl i tude 


340 


•*'40 


Amplitude 


10 ns - 1.09 ms 


KOR 1030 


3.0 -3.5 
(500) 


1, 

Phase 


182 


■*•10 


Doppler, 

Frequency 


<95 ns - 2.9 ms 


WHITTAKER 
MIP 830 


6 - 18 
(220) 


8, 

Amplitude 


1048 


-^3 


*** 


10 ns - 2.048 ms 
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V CONCLUSIONS AND RECOMMENDATIONS 

A. CONCLUSIONS 

This thesis presents a method to calculate the level of 
harmonics generated in the DRFM due to the sampling and 
quantization process. Fourier series analysis is used for the 
calculation of the harmonic levels. The quantization process 
will generate the harmonics of the fundamental signal. The 
sampling process will fold these harmonics into the baseband 
(IBW) of the DRFM as inter-modulation products of the sampling 
frequency (clock frequency of the A/D) and the signal 
frequency. The level of the harmonics is dependent on the 
number of quantization bits of the analog-to-digital 
converter. Further, it is also dependent on the ratio of the 
input signal frequency and the sampling frequency (normalized 
signal frequency) . If the normalized frequency is low, the 
study indicated that the harmonic levels in the output are low 
and correspondingly, if the normalized frequency is high, 
the harmonic level is high. 

The results obtained here are more conservative since the 
model used does not consider circuit limitations like 

• circuit stray capacitances which round off the quantized 
waveform 

• generation of even harmonics due to variations in 
components values, thus reducing odd harmonics 
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harmonic power spread due to local oscillator drifts 



The reduction of harmonic levels in single bit and multi 
bit DRFMs is possible. This can be accomplished by 
introducing time jitter in the sampling signal. Another method 
is to vary the duty cycle of the quantized output waveform. 

The harmonic computation programs written in MATLAB can 
accept number of bits as input and plot the harmonic power 
level vs harmonic number to a maximum of 99 harmonics. The 
total harmonic power is also computed by the program. The 
sampled and quantized signal harmonic power computation 
programs accept normalized input frequency with respect to 
sampling frequency. The plot and data for each harmonic level 
is generated. The user of these programs can determine spur 
levels of DRFM having any number of bits. 

In system applications, DRFMs offer many advantages over 
FMLs as deception jamming devices. The DRFM technology can be 
utilized to build broad bandwidth frequency sources for radar 
and communication applications. Also this technology can be 
used to store signals which can be analyzed either in real 
time or off-line. The DRFMs have been used in radar 
simulators for generation of different types of signals like 
target return, clutter signal, and jamming signal. 
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B . RE COMMENDAT IONS 



The DRFM is a new area of interest to the EW community. 
The reduction in harmonic levels due to sampling time jitter 
and duty cycle modulation should be investigated further. 
Phase encoded DRFMs and their characteristics have to be 
investigated as they offer some advantages over amplitude 
encoded DRFMs. Studies can also be conducted in the area of 
output frequency accuracy with respect to stored pulse length 
and sampling rate. The digital methods for imparting 
modulation to the acquired signal is another area of interest. 
Investigations are also suggested in the time multiplexed and 
frequency multiplexed techniques in DRFM. 
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APPENDIX A 



PLOTS OF MULTIBIT QUANTIZED SIGNAL 



The Appendix A contains the MATLAB program to calculate 
the harmonic of the multi bit quantized signal and the plots 
of harmonic power vs harmonic number for one to six bits 
quantized signal. The comparison of the harmonic power level 
of one to six bit amplitude quantized signal is given in 
Table A. 7. 
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ONE BIT AMPLITUDE QUANTIZED SIGNAL 
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Figure A.l One bit amplitude quantized signal 
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Figure A. 2 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of one bit amplitude quantized 
signal (expanded view, harmonics 1 to 19) 
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ONE BIT AMPLITUDE QUANTIZED SIGNAL 
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Figure A. 3 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of one bit amplitude quantized 
signal (harmonics 1 to 99) 



TABLE A.l HARMONIC POWER (IN DB) RELATIVE TO FUNDAMENTAL VS 
HARMONIC NUMBER OF ONE BIT AMPLITUDE QUANTIZED SIGNAL 



Harmonic 

Number 


Power 

level 


Harmonic 

number 


Power 

level 


Harmonic 

number 


Power j 
level 


1 


0.0 


19 


- 22.28 


29 


- 27.96 ' 


3 


- 9.54 


15 


- 23.52 


27 


- 28.63 


5 


- 13.98 


17 


- 24.61 


29 


- 29.25 


7 


- 16.90 


19 


- 25.58 


31 


- 29.83 


9 


- 19.08 


21 


- 26.44 


33 


- 30.37 


11 


- 20.83 


23 


- 27.23 


35 


- 30.88 
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TWO BIT AMPLITUDE QUANTIZED SIGNAL 
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Figure A. 5 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of two bit amplitude quantized 
signal (expanded view, harmonics 1 to 19) 
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Figure A. 6 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of two bit amplitude quantized 
signal (harmonics 1 to 99) 



TABLE A. 2 HARMONIC POWER (IN DB) RELATIVE TO FUNDAMENTAL VS 
HARMONIC NUMBER OF TWO BIT AMPLITUDE QUANTIZED SIGNAL 



Harmonic 

Number 


Power 

level 

(dB) 


Harmonic 

number 


Power 

level 

(dB) 


Harmonic 

number 


Power 

level 

(dB) 


1 


0 


13 


- 22.28 


25 


- 27.96 


3 


- 14.96 


15 


- 28.94 


27 


- 34.05 


5 


- 36.86 


17 


- 47.49 


29 


- 52.13 


7 


- 39.78 


19 


- 48.45 


31 


- 52.71 


9 


- 24.50 


21 


- 31.86 


33 


- 35.79 


11 


- 20.83 


23 


- 27.23 


35 


- 30.88 
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THREE BIT AMPLITUDE QUANTIZED SIGNAL 



Reconstructed signal 3 bits 




Figure A. 7 Three bit amplitude quantized signal 
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Figure A. 8 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of three bit amplitude 
quantized signal (expanded view, harmonics 1 to 19) 
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THREE BIT AMPLITUDE QUANTIZED SIGNAL 
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Figure A. 9 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of three bit amplitude 
quantized signal (harmonics 1 to 99) 



TABLE A. 3 HARMONIC POWER (IN DB) RELATIVE TO FUNDAMENTAL VS 
HARMONIC NUMBER OF THREE BIT AMPLITUDE QUANTIZED SIGNAL 



1 

Harmonic 

Number 


Power 

level 

(dB) 


Harmonic 

number 


Power 

level 

(dB) 


Harmonic 

number 


Power 

level 

(dB) 


! 1 


0.0 


13 


- 34.01 


25 


- 32.14 


3 


- 21.33 


15 


- 32.87 


27 


- 37.62 


5 


- 59.65 


17 


- 40.72 


29 


- 36.63 


7 


- 28.90 


15 


- 39.59 


31 


- 45.78 


9 


- 34.79 


21 


- 31.01 


33 


- 48.68 


11 


- 55.51 


23 


- 27.10 


35 


- 42.98 
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FOUR BIT AMPLITUDE QUANTIZED SIGNAL 
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Figure A. 10 Four bit amplitude quantized signal 
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Figure A. 11 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of four bit amplitude quantized 
signal (expanded view, harmonics 1 to 19) 



64 



FOUR BIT AMPLITUDE QUANTIZED SIGNAL 



Harmonic power Vs Harmonic number 
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Figure A. 12 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of four bit amplitude quantized 
signal (harmonics 1 to 99) 



TABLE A. 4 HARMONIC POWER (IN DB) RELATIVE TO FUNDAMENTAL VS 
HARMONIC NUMBER OF FOUR BIT AMPLITUDE QUANTIZED SIGNAL 



Harmonic 

Number 


Power 

level 

(dB) 


Harmonic 

number 


Power 

level 

(dB) 


Harmonic 

number 


Power 1 
level i 
(dB) 


1 


0.0 


19 


- 43.02 


25 


- 61.30 1 


3 


- 28.00 


19 


- 70.51 


27 


- 53.63 


5 


- 56.18 


17 


- 40.32 


25 


- 38.21 


7 


- 33.26 


19 


- 54.32 


31 


- 67.44 


9 


- 62.79 


21 


- 84.92 


33 


- 41.96 


11 


- 39.96 


23 


- 39.79 


35 


- 40.67 



65 



FIVE BIT AMPLITUDE QUANTIZED SIGNAL 



Reconstruct«scJ signal 5 bits 




Figure A. 13 Five amplitude reconstructed signal 
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Figure A. 14 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of five bit amplitude quantized 
signal (expanded view, harmonics 1 to 19) 
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FIVE BIT AMPLITUDE QUANTIZED SIGNAL 
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Figure A. 15 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of five bit amplitude quantized 
signal, (harmonics 1 to 99) 



TABLE A. 5 HARMONIC POWER (IN DB) RELATIVE TO FUNDAMENTAL VS 
HARMONIC NUMBER OF FIVE BIT AMPLITUDE QUANTIZED SIGNAL 



Harmonic 

Number 


Power 

level 

(dB) 


Harmonic 

number 


Power 

level 

(dB) 


Harmonic 

number 


Power 

level 

(dB) 


1 


0.00 


13 


- 69.32 


25 


- 49.24 


3 


- 34.72 


15 


- 47.79 


27 


- 62.70 


5 


- 52.35 


17 


- 54.30 


29 


- 63.65 


7 


- 39.67 


19 


- 58.01 


31 


- 49.97 


9 


- 86.66 


21 


- 49.28 


33 


- 57.84 


11 


- 43.02 


23 


- 71.22 


35 


- 50.92 
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SIX BIT AMPLITUDE QUANTIZED SIGNAL 



Reconstructed sf^nol 6 bits 




Figure A. 16 Six bit amplitude quantized signal 
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Figure A. 17 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of six bit amplitude quantized 
signal (expanded view, harmonics 1 to 19) 
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SIX BIT AMPLITUDE QUANTIZED SIGNAL 



Harmonic power Vs Harmonic number 
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Figure A. 18 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of six bit amplitude quantized 
signal (harmonics 1 to 99) 



TABLE A. 6 . HARMONIC POWER (IN DB) RELATIVE TO FUNDAMENTAL VS 
HARMONIC NUMBER OF SIX BIT AMPLITUDE QUANTIZED SIGNAL 



Harmonic 

Number 


Power 

level 

(dB) 


Harmonic 

number 


Power 

level 

(dB) 


Harmonic 

number 


Power 

level 

(dB) 


1 


0.00 


13 


- 92.20 


25 


- 61.96 


3 


- 41.39 


15 


- 51.56 


27 


- 67.19 


5 


- 54.23 


17 


- 89.04 


29 


- 58.55 


7 


- 46.57 


IS 


- 54.49 


31 


- 93.48 


9 


- 68.74 


21 


- 69.00 


33 


- 57.82 


11 


- 49.22 


23 


- 58.95 


35 


- 84.77 
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TABLE A. 7 HARMONIC POWER (DB) VS HARMONIC NUMBER 



Harmonic 

Number 


— 

Power level in dB of the Harmonics relative to 
the fundamental (Quantization bits) 


1 bit 


2 bit 


3 bit 


4 bit 


5 bit 


6 bit 


1 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


3 


- 9,54 


- 14.96 


- 21.33 


- 28.00 


- 34.72 


- 41.39 


5 


- 13.98 


- 36.86 


- 59.65 


- 56.18 


- 52.35 


- 54.23 


7 


- 16.90 


- 39.78 


- 28.90 


- 33.26 


- 39.67 


- 46.57 


9 


- 19.08 


- 24.50 


- 34.79 


- 62.70 


- 86.66 


- 68.74 


11 


- 29.83 


- 20.83 


- 55.51 


- 39.96 


- 43,02 


- 49.22 


13 


- 22.28 


- 22.28 


-34.01 


- 43.02 


- 69.32 


- 92.20 


15 


- 23.52 


- 28.94 


- 32.87 


- 70.51 


- 47.79 


- 51.56 


17 


- 26.44 


- 47.49 


- 40.72 


- 40.32 


- 51.30 


- 89.04 


19 


- 25.58 


- 48.45 


- 39.96 


- 54.32 


- 58.01 


- 54.49 


21 


- 26.44 


- 31.86 


- 31.01 


- 84.92 


- 49.28 


- 69.00 


23 


- 27.23 


- 27.23 


- 27.10 


- 39.79 


- 71.22 


- 58.95 


25 


- 27.96 


- 27.96 


- 32.14 


- 61.30 


- 49,94 


- 61.96 


27 


- 29.83 


- 34.05 


- 37.62 


- 53.63 


- 62.70 


- 67.19 


25 


- 29.25 


- 52.13 


- 36.63 


- 38.21 


- 63.65 


- 58.55 


31 


- 29.83 


- 52.71 


- 45.78 


- 67.44 


- 49.97 


- 93.48 


33 


- 30.37 


- 35.79 


- 48.68 


- 41.96 


- 57.84 


- 57.82 


29 


- 30.88 


- 30.88 


- 42.98 


- 40.67 


- 50.92 


- 84.77 
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A. MATLAB PROGRAM FOR CALCULATION OF HARMONIC POWER IN MULTI 
BIT QUANTIZATION. 

% FILE NAME 
% OPERATING SYSTEM 
% SYSTEM 
% 

% M BIT NOOFLEVELS 

!del ampll.met 
!del ampl2.met 
!del ampl3.met 

!del ampll.pic 
!del ampl2.pic 
!del amplS.pic 

clear 

axis ( ' square ' ) 

bits = input ('Enter number of bits 1 to 8 = '); % no of bits 
N = (2^bits)/2; % no of levels 

AO = 1; % amplitude of the sinusoid 

Maxhar =50; % calculation is done up to maximum 

% of 99 harmonics; 

A2 = 4*A0/(pi*N); 

% HARMONIC CONTENT CALCULATION 

for k = 1;N; % STEP NO IN HALF CYCLE 

for n = 1; Maxhar, 

s(n) = cos ( (2*(n-l)+l) * asin ( (k-1) /(N) ) ); 
hno(n) = 2*(n-l)+l; 
fl(n,k) = A2/(2*(n-l)+l) *s(n) ; 
theta(n) = (2* (n-1) +1) * (2*pi) ; 
end, 
end, 

disp (' Harmonic calculation complete') 

pnk=zeros (1 :hno) ; 
for k = 1:N 

pnk = pnk + fl(:,k); 

end. 



; DRQUA.M 
; DOS 5 . 0 
: PC 486/33 

= (2^M)/2 WHERE M IS NO OF BITS 
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disp (' calculating the log') 
pk = pnk . * pnk ; 
p = pk / max(pk) ; 
pp = 10*logl0(p); 



% HARMONIC POWER CALCULATIONS 
disp (' calculating max harmonic') 
horpow = 0; 
for i = 2:Maxhar, 

horpow = horpow + 10^ (pp (i)/10) ; 
end, 

horpow4bit = 10*logl0 (horpow) ; 



textl = 
text2 = 
texts = 
texts = 
text4 = 
text4 = 
texts = 
texts = 
texts = 
text? = 
texts = 
texts 1 = 



[ 'kf21/nityan/dramp.m' ] ; 

[num2str (bits) ' bit Amplitude quantization']; 
['Total harmonic power : ']; 

[texts num2str (horpow4bit) 'dBr']; 

['Maximum harmonic power at ']; 

[text4 num2str (hno(2) ) 'rd harmonic']; 

[num2str (hno(2) ) ' rd harmonic power level : ']; 

[texts num2str (pp(2) ) 'dBr']; 

['Harmonic power Vs Harmonic number']; 

['Relative to fundamental < ']; 

[ 'Harmonic number' ] ; 

[ 'Normalized frequency' ] ; 



% TIME PLOT OF THE RECONSTRUCTED SIGNAL 

disp (' calculating reconstructed signal ') t=0:pi/2S6: 2*pi 
for z= 1 ; length (t), 

fz= pnk .* sin (theta' * t(z)); 
v=0 ; 

for u= 1 : length (fz) 
v=v+f z (u) ; 
end, 

ss(z) = v; 
end, 

axis([0 7 -l.S*max(ss) l.S*max(ss) ]) ; 
plot (t , ss) 
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if bits == 1 

title ([ 'Reconstructed signal ' num2str (bits) 

else 

title ([ 'Reconstructed signal ' num2str (bits) 

end 

xlabel ( ' Time ' ) 
ylabel ( 'Amplitude ' ) 
meta ampll 



% PLOT THE DATA 
axis 

plot(hno(l: 10) ,pp(l: 10) ) 

title (text6) 

ylabel (text?) 

xlabel (texts) 

grid 



text (0.4, 0.85, 
text (0 . 4,0.80, 
text (0.4,0.75, 
text (0.4,0.70, 
meta amp 12 

plot (hno,pp) 
title (text6) 
ylabel (text7) 
xlabel (texts) 
grid 

text (0.4,0.85, 
text (0.4,0.80, 
text (0.4,0.75, 
text (0.4, 0.70, 
meta amp 13 



text2, 'sc') 
text4 , 'sc') 
text5, 'sc') 
texts, 'sc') 



text2 , 'sc') 
text4, 'sc') 
text5, 'sc') 
texts, 'sc') 



disp ( 'Writing data to the file QUA.dat') 
filename = 'qua.dat' ; 

noofbits = [' No of bits = ' num2str (bits) '\n' 
fprintf (filename, noofbits); 

fprintf( filename, ' Harmonic Power \n' 

f printf ( filename, ' number level \n ' 



' bit']) 

' bits ' ] ) 



73 



'%6.0f %12.2f \n', hno(k), pp(k)) 



for k = 1:20, 

fprintf (filename , 
end, 

!gpp ampll /dhpgl /fampll.pic 
!gpp ampl2 /dhpgl /fampl2.pic 
!gpp ampl3 /dhpgl /fampl3.pic 
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APPENDIX B PLOTS OF ERROR MINIMIZED QUANTIZED SIGNAL 



Appendix B contains the MATLAB programs to calculate the 
harmonic power of error minimized 2, 3 ,4 and 5 bit amplitude 
quantized signal. The plots of harmonic power vs harmonic 
number for error minimized quantized signal are also presented 
in this appendix. 
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TWO BIT AMPLITUDE QUANTIZED SIGNAL (ERROR MINIMIZED) 



pow«r vo harmonlo Number 




Figure B.l Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of two bit amplitude quantized 
signal (squared error minimized, expanded view, harmonics 1 to 
19) 



Harmonic power vs Hormonlo Number 




Harmonic Number 



Figure B.2 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of two bit amplitude quantized 
signal (squared error minimized, harmonics 1 to 99) 
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THREE BIT AMPLITUDE QUANTIZED SIGNAL (ERROR MINIMIZED) 



<1_> 




Hormonl^ pow^r vo hormonfo Numt>«r 




Hormonic numtxsr 



Figure B.3 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of three bit amplitude 
quantized signal (squared error minimized, expanded view, 
harmonics 1 to 19) 



Harmonic power vs Harmonic Numt>er 




Figure B.4 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of three bit amplitude 
quantized signal (squared error minimized, harmonics 1 to 99) 
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FOUR BIT AMPLITUDE QUANTIZED SIGNAL (ERROR MINIMIZED) 



H^rmonlo pow^r vo hormpnlp Numfc>«r 




Figure B.5 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of four bit amplitude quantized 
signal (squared error minimized, expanded view, harmonics 1 to 
19) 



Harmonic power vs harmonic Number 




Figure B.6 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of four bit amplitude quantized 
signal (squared error minimized, harmonics 1 to 99) 
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FIVE BIT AMPLITUDE QUANTIZED SIGNAL (ERROR MINIMIZED) 



pow^r \/o hormprifo Mumb-cr- 




Figure B.7 Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of five bit amplitude quantized 
signal (squared error minimized, expanded view, harmonics 1 to 
19) 



Horrmonlc power vs hormonlo Number 




nurnb^r 



Figure B.8. Plot of harmonic power (in dB) relative to 
fundamental vs harmonic number of five bit amplitude quantized 
signal (squared error minimized, harmonics 1 to 99 plotted) 
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A. MATLAB PROGRAM FOR CALCULATION OF HARMONIC POWER IN 2 BIT 



AMPLITUDE QUANTIZED (ERROR MINIMIZED) SIGNAL 



% 



FILE NAME 
OPERATING SYSTEM 
SYSTEM 



MODBIT2.M 
DOS 5.0 
PC 486/33 



% 



% 



clear 

!del modbit21 .met 
!del modbit22.met 

% SET THE AMPLITUDE STEP VALUE AND TRANSITION ANGLES 
a = 0.4188161; 
t2 = asin(3*a/2); 

% TWO BIT HARMONIC CALCULATION 

i=0: 100; 

n=2*i+l; 

for i = l;length(n), 
a2 = cos(n(i) *t2) ; 
s(i) = 4*a/ (n(i) *pi) * (l+a2) ; 

end 

p = 10*logl0(s .* s) - 10*logl0 (s (1) * s(l)); 



% CALCULATION OF TOTAL HARMONIC POWER 

harpower = 0 ; 

for i = 2: length (p) 

harpower = harpower +10^ (p (i) /lO) ; 

end, 

harpow = harpower*100 ; 

% PLOT THE DATA 

ttext = 'Harmonic power vs harmonic Number' ; 
xtext = 'Harmonic Number'; ytext = 'Harmonic Power'; 
text2 = 'Two bit amplitude quantization'; 
text! = 'Maximum harmonic power at 9th harmonic'; 
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text4 = ['9th Harmonic Power level = num2str (p(5) ) , 'dBr ' ] 
texts = ['Total harmonic power ' , num2str (harpow) , '%']; 

plot (n(l: 10) , p(l:10)); title (ttext) ; 
xlabel (xtext) ; ylabel (ytext) ; grid 
text (0 . 3 , 0 . 85 , text2 , 'sc'); 

text(0.3,0.80, text3, 'sc'); 
text (0.3, 0.75, text4, 'sc'); 
text (0.3, 0.70, texts, 'sc') 
meta modbit21; pause; 

!gpp modbit21 /dhpgl /fmodbit21.pic 

plot (n ( 1 : 50) , p(l:50)); title (ttext) ; 
xlabel (xtext) ; ylabel (ytext) ; grid 

text (0.3, 0.85, text2,'sc'); text (0 . 3 , 0 . 80 , text 3, 'sc') 
text (0 . 3 , 0 . 75, text4 , 'sc'); text (0 . 3 , 0 . 70 , texts, 'sc') 

meta modbit22; pause; 

!gpp modbit22 /dhpgl /fmodbit22 . pic 
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B. MATLAB PROGRAM FOR CALCULATION OF HARMONIC POWER IN 3 BIT 
AMPLITUDE QUANTIZED (ERROR MINIMIZED) SIGNAL 

% FILE NAME : MODBIT3.M 

% OPERATING SYSTEM : DOS 5.0 

% SYSTEM : PC 486/33 

clear 

!del modbit31 .met 
!del modbit32.met 

% AMPLITUDE STEP a 
a = 0.2332685; 
t2 = asin(3*a/2) ; 
t3 = asin(5*a/2) ; 
t4 = asin(7*a/2) ; 

% THREE BIT HARMONIC CALCULATIONS 
i=0: 100; 
n=2*i+l ; 

for i = l:length(n), 
a2 = cos(n(i) *t2) ; 
a3 = cos (n (i) *t3 ) ; 
a4 = cos(n(i) *t4) ; 

s(i) = 4*a/ (n (i) *pi) * (l+a2+a3+a4) ; 

end 

p = 10*logl0(s .* s) - 10*logl0(s(l) * s(l)); 

% CALCULATION OF THE TOTAL HARMONIC POWER 

harpower = 0 ; 

for i = 2: length (p) 

harpower = harpower +10^ (p (i)/10) ; 

end, 

harpow = harpower* 10 0 ; 
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% PLOT THE DATA 



ttext = 
xtext = 
ytext = 
textl = 
text2 = 
texts = 
texts = 
text4 = 



•Harmonic power vs harmonic Number'; 

•Harmonic number'; 

•Harmonic power'; 

•Three bit amplitude quantization'; 

•Maximum harmonic power at 19th harmonic'; 

' 19th Harmonic Power level = ' ; 

[texts numSstr (p ( 10) ) 'dBr']; 

['Total harmonic power ' , numSstr (harpow) , '%']; 



plot (n(l: 10) , p(l;10)) 

title (ttext) ; xlabel (xtext) ; ylabel (ytext) ; grid 

text (0 . S , 0 . 85 , textl, 'sc'); 

text (0. S , 0 . 80 , texts, 'sc') 

text (0. S , 0 . 75 , texts, 'sc'); 

text (0. S , 0 . 70 , text4 , 'sc') 

meta modbitSl; pause 

plot(n(l:50) , p(l:50)) 

title (ttext) ; xlabel (xtext) ; ylabel (ytext) ; grid 

text (0 . S , 0 . 85 , textl, 'sc') 

text (O.S, 0.80, texts, 'sc') 

text (0 . S , 0 . 75, texts, 'sc') 

text (0 • S , 0 . 70 , text4 , 'sc') 

meta modbitSS 

!gpp modbitSl /dhpgl /fmodbitSl .pic 
!gpp modbitSS /dhpgl /fmodbitSS .pic 
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C. MATLAB PROGRAM FOR CALCULATION OF HARMONIC POWER IN 4 BIT 

AMPLITUDE QUANTIZED (ERROR MINIMIZED) SIGNAL 

% FILE NAME : MODBIT4.M 

% OPERATING SYSTEM : DOS 5.0 

% SYSTEM : PC486/33 

clear 

!del modbit41.met 
!del modbit42 .met 

% SET AMPLITUDE STEP VALUE AND TRANSITION ANGLES 

a = 0.12139; 

t2 = asin(3*a/2) ; 

t3 = asin(5*a/2); 

t4 = asin(7*a/2); 

t5 = asin(9*a/2); 

t6 = asin(ll*a/2) ; 

t7 = asin(13*a/2) ; 

t8 = asin(15*a/2) ; 



% THREE BIT HARMONICS CALCULATIONS 
i=0: 100; 



n=2*i+l 
for i = 
a2 = 
a3 = 
a4 = 
a5 = 
a6 = 
a7 = 
a8 = 
s(i) 

end 



1 : length (n) , 
cos(n(i) *t2) ; 
cos(n(i) *t3) ; 
cos (n (i) *t4 ) ; 
cos(n(i) *t5) ; 
cos(n(i) *t6) ; 
cos(n(i) *t7) ; 
cos (n(i) *t8) ; 

= 4*a/ (n(i) *pi) 



* ( l+a2+a3+a4+a5+a6+a7+a8 ) 



p = 10*logl0(s .* s) - 10*logl0 (s (1) * s(l)); 
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% CALCULATION OF TOTAL HARMONIC POWER 

harpower = 0; 

for i = 2: length (p) 

harpower = harpower +10^ (p (i) /lO) ; 

end, 

harpow = harpower*100 ; 

% PLOT THE DATA 

ttext = 'Harmonic power vs harmonic Number'; 

xtext = 'Harmonic number'; 

ytext = 'Harmonic power'; 

textl = 'Four bit amplitude quantization'; 

text2 = 'Maximum harmonic power at 43rd harmonic'; 

texts = ['43th Harmonic Power level = ' num2str (p(22) ) 'dBr'] 

text4 = ['Total harmonic power ', num2str (harpow) , '%']; 



plot (n(l: 10) , p(l:10)) 

title (ttext) ; xlabel (xtext) ; ylabel (ytext) ; grid 
text (0. 3 , 0 . 85, textl, 'sc'); text (0. 3 , 0 . 80, text2 , 'sc') 

text (0. 3, 0.75, texts, 'sc'); text ( 0 . 3 , 0 . 70 , text4 , 'sc') 

meta modbit41; pause 

plot(n(l:50) , p(l:50)) 
title (ttext) ; xlabel (xtext) ; 
text ( 0 . 3 , 0. 85, textl, 'sc'); 
text ( 0 . 3 , 0 . 75 , text3 , 'sc'); 

meta modbit42 

Igpp modbit41 /dhpgl /fmodbit41.pic 
!gpp modbit42 /dhpgl /fmodbit42 .pic 



ylabel (ytext) ; grid 

text (0.3,0.80, texts , ' sc ' ) 

text ( 0. 3 , 0 . 70, text4, 'sc') 
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D. MATLAB PROGRAM FOR CALCULATION OF HARMONIC POWER IN 5 BIT 
AMPLITUDE QUANTIZED (ERROR MINIMIZED) SIGNAL 



!del modbitSl .met 
!del modbit52.met 

clear 

%SET THE AMPLITUDE STEP AND TRANSITION ANGLE VALUES 

a = 0.0616995; 

t2 = asin(3*a/2); 

t3 = asin(5*a/2); 

t4 = asin(7*a/2) ; 

t5 = asin(9*a/2); 

t6 = asin(ll*a/2) ; 

t7 = asin(13*a/2) ; 

t8 = asin(15*a/2) ; 

t9 = asin(17*a/2) ; 

tlO = asin(19*a/2) ; 

til = asin(21*a/2) ; 

tl2 = asin(23*a/2) ; 

tl3 = asin(25*a/2) ; 

tl4 = asin(27*a/2) ; 

tl5 = asin(29*a/2) ; 

tl6 = asin(31*a/2) ; 



% FILE NAME 
% OPERATING SYSTEM 
% SYSTEM 



MODBIT5.M 
DOS 5.0 
PC486/33 
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% FIVE BIT BIT HARMONIC CALCULATIONS 

i=0:100; 

n=2*i+l ; 

for i = l:length(n), 
a2 = cos (n ( i) *t2 ) ; 
a3 = cos (n (i) *t3 ) ; 
a4 = cos (n (i) *t4 ) ; 
a5 = cos(n(i) *t5) ; 
a6 = cos (n (i) *t6) ; 
a7 = cos (n (i) *t7 ) ; 
a8 = cos (n (i) *t8 ) ; 
a9 = cos(n(i) *t9) ; 
alO = cos (n(i) *tlO) ; 
all = cos(n(i) *tll) ; 
al2 = cos (n (i) *tl2 ) ; 
al3 = cos (n (i) *tl3 ) ; 
al4 = cos (n (i) *tl4) ; 
al5 = cos (n (i) *tl5) ; 
al6 = cos (n (i) *tl6) ; 

a = l+a2+a3+a4+a5+a6+a7+a8+a9 # 
a = a+al0+all+al2+al3+al4+al5+al6 ; 
s(i) = 4*a/(n(i)*pi)*(a) ; 

end 

p = 10*logl0(s .* s) - 10*logl0 (s(l) * s(l)); 



% CALCULATION OF TOTAL HARMONIC POWER 

harpower = 0 ; 

for i = 2: length (p) 

harpower = harpower +10^ (p (i)/10) ; 

end, 

harpow = harpower*? 00 ; 
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% PLOT THE DATA 



ttext 

xtext 

ytext 

textl 

text2 

texts 

texts 

text4 



'Harmonic power vs harmonic Number'; 
'Harmonic number' ; 

'Harmonic power' ; 

'Five bit amplitude quantization'; 

'Maximum harmonic power at 80th harmonic'; 
['80th Harmonic Power level = ']; 

[texts numSstr (p(40) ) 'dBr']; 

['Total harmonic power ' , numSstr (harpow) , 



plot(n(l:10) , p(l:10)) 

title (ttext) ; xlabel (xtext) ; ylabel (ytext) ; grid 
text (O.S, 0.85, textl, 'sc'); text (0 , S , 0 . 80 , texts, 
text ( 0 . S , 0 . 75 , texts, 'sc'); text (0 . S , 0 . 70 , text4 , 
meta modbit51; pause 



plot (n(l: 50) , p(l:50)) 

title (ttext) ; xlabel (xtext) ; ylabel (ytext) ; grid 
text (0 . S , 0 . 85 , textl, ' sc ' ) ; text (0 . S , 0 . 80 , text2 , 
text ( 0. S, 0.75, texts, ' sc ' ) ; text (0 . S , 0 . 70 , text4 , 
meta modbit52 ; 



!gpp modbit51 /dhpgl 
!gpp modbit52 /dhpgl 



/fmodbit51 . pic 
/fmodbit52 . pic 



'sc' ) 
' sc ' ) 



'sc' ) 
' sc ' ) 
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APPENDIX C 



PLOTS OF SAMPLED AND QUANTIZED SIGNAL 



Appendix C contains all the plots of harmonic power vs 
normalized signal frequency for one to four bit of quantized 
and sampled signal. The MATLAB program for calculating the 
harmonic power level is also included in this appendix. 
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ONE BIT AMPLITUDE QUANTIZED SAMPLED SIGNAL 



1 bit omplitLJcJe cquantized Sc SBamplecJ signal 
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Figure C.l Plot of harmonic power (in dB) relative to 
fundamental vs normalized frequency of one bit amplitude 
quantized sampled signal. Normalized input signal 

frequency = 0.005. (Plot from dc to half the sampling 
frequency) 



T bit amplitude quantized Sc sampled signal 
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Figure C.2 Plot of harmonic power (in dB) relative to 
fundamental vs normalized frequency of one bit amplitude 
quantized sampled signal. Normalized input signal 

frequency = 0.025. (Plot from dc to half the sampling 
frequency) 
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ONE BIT AMPLITUDE QUANTIZED SAMPLED SIGNAL 



1 fc>it omplltucde: quantiz«<d £c samplecJ signal 




Figure C.3 Plot of harmonic power (in dB) relative to 
fundamental vs normalized frequency of one bit amplitude 
quantized sampled signal. Normalized input signal 

frequency = 0.05. (Plot from dc to half the sampling 
frequency) 
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Figure C.4 Plot of harmonic power (in dB) relative to 
fundamental vs normalized frequency of one bit amplitude 
quantized sampled signal. Normalized input signal frequency = 
0.1. (Plot from dc to half the sampling frequency) 
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TABLE C.l HARMONIC POWER W.R.T FUNDAMENTAL VS HARMONIC NUMBER 
OF ONE BIT SAMPLED AND QUANTIZED SIGNALS 



Harmonic 

Number 


Power level relative to fundamental in dB 


Sampled and quantized signal of 
normalized frequency 


Quantized 
signal only 
( frequency 
independent) 


Freq 
= 0.005 


Freq 
= 0.025 


Freq 
= 0.05 


Freq 
= 0.1 


1 


0.00 


0.00 


0.00 


0.00 


0.00 


3 


- 9.45 


- 8.86 


- 8.23 


- 7.01 


- 9.54 


5 


- 13.80 


- 12.64 


- 11.41 


***** 


- 13.98 


7 


- 16.64 


- 14.93 


- 13.05 


***** 


- 16.90 


9 


- 18.73 


- 16.49 


- 13.83 


***** 


- 19.08 


11 


- 20.38 


- 17.61 


***** 


***** 


- 20.83 


13 


- 21.74 


- 18.42 


***** 


***** 


- 22.28 


15 


- 22.89 


- 19.00 


***** 


***** 


- 23.52 


17 


- 23.88 


- 19.39 


***** 


***** 


- 24.61 


19 


' - 24.75 


- 19.61 


***** 


***** 


-25.S8 1 



NOTE; ***** INDICATES THAT THE HARMONIC IS 
BEYOND HALF THE SAMPLING FREQUENCY AND 
HENCE NEGLECTED. 
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TWO BIT AMPLITUDE QUANTIZED SAMPLED SIGNAL 
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Figure C.5 Plot of Harmonic power (in dB) relative to 
fundamental vs normalized frequency of two bit amplitude 
quantized sampled signal. Normalized input signal 

frequency = 0.005. (Plot from dc to half the sampling 
frequency) 
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Figure C.6 Plot of harmonic power (in dB) relative to 
fundamental vs normalized frequency of two bit amplitude 
quantized sampled signal. Normalized input signal frequency = 
0.025. (Plot from dc to half the sampling frequency) 
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TWO BIT AMPLITUDE QUANTIZED SAMPLED SIGNAL 
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Figure C.7 Plot of harmonic power (in dB) relative to 
fundamental vs normalized frequency of two bit amplitude 
quantized sampled signal. Normalized input signal frequency = 
0.05. (Plot from dc to half the sampling frequency) 
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Figure C.8 Plot of harmonic power (in dB) relative to 
fundamental vs normalized frequency of two bit amplitude 
quantized sampled signal. Normalized input signal 
frequency = 0.1. (Plot from dc to half the sampling frequency) 
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TABLE C.2 HARMONIC POWER W.R.T FUNDAMENTAL VS HARMONIC NUMBER 
OF TWO BIT SAMPLED AND QUANTIZED SIGNALS 



Harmonic 

Number 


Power level relative to fundamental in dB 


Sampled and quantized signal of 
normalized frequency 


Quantized 
signal only 
(frequency j 
independent) i 


Freq 
= 0.005 


Freq 
= 0.025 


Freq 
= 0.05 


Freq 
= 0.1 


1 


0.00 


0.00 


0.00 


0.00 


0.00 


9 


- 14.76 


- 13.44 


- 12.18 


- 9.94 


- 14.96 


5 


- 35.38 


- 25. 61 


- 21.55 


kkkkk 


- 36.86 


7 


- 36.27 


- 26.82 


- 20.41 


kkkkk 


- 39.78 


9 


- 23.81 


- 19.99 


- 16.13 


kkkkk 


- 24.50 


11 


- 20.60 


- 17.79 


•k-k-k-k-k 


kkkkk 


- 20.83 


13 


- 21.74 


- 19.05 


kk-kk-k 


kkkkk 


- 22.28 


15 


- 27.76 


- 21.69 


kkkkk 


kkkkk 


- 28.94 


17 


- 42.96 


- 25.61 


kkk kk 


kkkkk 


- 47.49 


" 1 


- 40.67 


- 26.89 


kkkkk 


kkkkk 


- 48.45 



NOTE: ***** INDICATES THAT THE HARMONIC IS 

BEYOND HALF THE SAMPLING FREQUENCY AND 
HENCE NEGLECTED. 
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THREE BIT AMPLITUDE QUANTIZED SAMPLED SIGNAL 
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Figure C.9 Plot of harmonic power (in dB) relative to 
fundamental vs normalized frequency of three bit amplitude 
quantized sampled signal. Normalized input signal 
frequency = 0.005. (Plot from dc to half the sampling 
frequency) 
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Figure C.IO Plot of harmonic power (in dB) relative to 
fundamental vs normalized frequency of three bit amplitude 
quantized sampled signal. Normalized input signal 
frequency = 0.025. (Plot from dc to half the sampling 
frequency) 
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THREE BIT AMPLITUDE QUANTIZED SAMPLED SIGNAL 
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Figure C.ll Plot of harmonic power (in dB) relative to 
fundamental vs normalized frequency of three bit amplitude 
quantized sampled signal. Normalized input signal 

frequency = 0.05. (Plot from dc to half the sampling 
frequency) 
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Figure C.12 Plot of harmonic power (in dB) relative to 
fundamental vs normalized frequency of three bit amplitude 
quantized sampled signal. Normalized input signal 
frequency = 0.1. (Plot from dc to half the sampling frequency) 
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TABLE C.3 HARMONIC POWER W.R.T FUNDAMENTAL VS HARMONIC NUMBER 
OF THREE BIT SAMPLED AND AMPLITUDE QUANTIZED SIGNALS 



Harmonic 

Number 


Power level relative to fundamental in dB 


Sampled and quantized signal of 
normalized frequency 


Quantized 
signal only 
( frequency 
independent) 


Freq 
= 0.005 


Freq 
= 0.025 


Freq 
= 0.05 


Freq 
= 0.1 


1 


0.00 


0.00 


0.00 


0.00 


0.00 


3 


- 20.86 


- 18.28 


- 16.34 


- 12.96 


- 21.33 


9 


- 54.70 


- 28.97 


- 22.23 


***** 


- 59.65 


1 


- 27.97 


- 22.95 


- 18.77 


***** 


- 28.90 


9 


- 33.23 


- 25.09 


- 18.25 


***** 


- 34.79 


11 


- 51.85 


- 28.42 


***** 


***** 


- 55.51 


13 


- 32.30 


- 24.68 


***** 


***** 


- 34.01 


15 


- 30.93 


- 23.23 


***** 


***** 


- 32.87 


17 


- 44.24 


- 27.85 


***** 


***** 


- 40.72 


19 


- 44.74 


- 30.18 


***** 


***** 


- 39.59 



NOTE: ***** INDICATES THAT THE HARMONIC IS 

BEYOND HALF THE SAMPLING FREQUENCY AND 
HENCE NEGLECTED. 
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FOUR BIT AMPLITUDE QUANTIZED SAMPLED SIGNAL 
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Figure C.13 Plot of harmonic power (in dB) relative to 
fundamental vs normalized frequency of four bit amplitude 
quantized sampled signal. Normalized input signal 

frequency = 0.005. (Plot from dc to half the sampling 
frequency) 
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Figure C.14 Plot of harmonic power (in dB) relative to 
fundamental vs normalized frequency of four bit amplitude 
quantized sampled signal. Normalized input signal 
frequency = 0.025. (Plot from dc to half the sampling 
frequency) 
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FOUR BIT AMPLITUDE QUANTIZED SAMPLED SIGNAL 
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Figure C.15 Plot of harmonic power (in dB) relative to 
fundamental vs normalized frequency of four bit amplitude 
quantized sampled signal. Normalized input signal 

frequency = 0.05. (Plot from dc to half the sampling 
frequency) 
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Figure C.16 Plot of harmonic power (in dB) relative to 
fundamental vs normalized frequency of four bit amplitude 
quantized sampled signal. Normalized input signal 

frequency = 0.1. (Plot from dc to half the sampling frequency) 
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TABLE C.4 HARMONIC POWER W.R.T FUNDAMENTAL VS HARMONIC NUMBER 
OF FOUR BIT SAMPLED AND QUANTIZED SIGNALS 



Harmonic 

Number 


Power level relative to fundamental in dB 


Sampled and quantized signal of 
normalized frequency 


Quantized 
signal only 
( frequency 
independent) 


Freq 
= 0.005 


Freq 
= 0.025 


Freq 
= 0.05 


Freq 
= 0.1 


1 


0.00 


0.00 


0.00 


0.00 


0.00 


3 


- 26.98 


- 22.66 


- 19.10 


- 15.23 


- 28.00 


5 


- 44.79 


- 28.81 


- 21.99 


***** 


- 56.18 


7 


- 32.12 


- 24.33 


- 20.35 


***** 


- 33.26 


9 


- 46.72 


- 28.54 


- 20.99 


***** 


- 62.79 


11 


- 37.03 


- 26.07 


***** 


***** 


- 39.96 


13 


- 39.12 


- 26. S8 


***** 


***** 


- 43.02 


15 


- 45.84 


- 26.88 


***** 


***** 


- 70.51 


17 


- 37.62 


- 24.35 


***** 


***** 


- 40.32 


19 


- 43.36 


- 27.22 


***** 


***** 


- 54.32 



NOTE: ***** INDICATES THAT THE HARMONIC IS 

BEYOND HALF THE SAMPLING FREQUENCY AND 
HENCE NEGLECTED. 
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A. MATLAB PROGRAM FOR CALCULATION OF HARMONIC POWER OF 
QUANTIZED SAMPLED SIGNAL. 

% FILENAME : DRAMP.M 

% OPERATING SYSTEM : DOS 5.0 
% SYSTEM : PC 486/33 

% M BIT NOOFLEVELS = (2^M)/2 WHERE M IS NO OF BITS 



clear 

!del drl.met 

bits = input (' No of bits of quantization (1 to 8) = '); 

N = (2^bits)/2; % no of levels 

AO = 1; 

fo= input ('Input signal frequency = '); 

To = 1/fo; 
wo = 2*pi*fo; 



Ac = 1; 
f c = 1 ; 

Tc = 1/fc; 
wc = 2*pi*fc; 
tc = 0.001*Tc; 

Maxhar = fix(fc/fo) ; 

harp = Maxhar/2 ; 

A1 = Ac*tc/Tc; 

A2 = 4*A0/(pi*N); 



% CALCULATION IS DONE UP TO MAXIMUM 
% HARMONICS; 



textl = [ 'kf21/nityan/dramp.in' ] ; 

text2 = [num2str (bits) 'bit Amplitude quantization']; 

% HARMONIC CONTENT CALCULATION 

for k = 1:N; % step no in half cycle 

for n = 1: Maxhar, 

s(n) = cos ( (2*(n-l)+i) * asin ( (k-1) /(N) ) ); 

hno(n) = 2*(n-l)+l; 

fl(n,k) = Al*A2/(2*(n-l)+l) *s(n) ; 
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theta(n) = (2* (n-1) +1) * (2*pi*fo) ; 
end, 
end, 

disp ( •Harmonic calculation complete') 



pnk=zeros (1 ;hno) ; 
for k = 1;N 

pnk = pnk + fl(:,k); 

end, 

disp (' calculating the log') 
pk = pnk .* pnk; 
p = pk / max(pk) ; 
pp = 10*logl0 (p) ; 

% HARMONIC POWER CALCULATIONS 
disp (' calculating max harmonic') 
harpow = 0 ; 
for i = 2:Maxhar, 

harpow = harpow + 10^ (pp (i)/10) ; 
end, 

harpow4bit = 10*logl0 (harpow) ; 



texts = 
texts = 
text4 = 
text4 = 
texts = 
texts = 
texts = 
textV = 
texts = 
texts 1 = 
textsp = 
textlO = 
textlO = 
textll = 
textll = 
textlS = 
textlS = 



['Total Harmonic power : ']; 

[texts num2str (harpow4bit) 'dBr']; 

['Max harmonic power at ']; 

[text4 num2str (hno(2) ) 'rd harmonic']; 

[num2str (hno (2) ) 'rd Harmonic power level : ']; 

[texts num2str (pp(2) ) 'dBr']; 

['Harmonic power Vs Harmonic number']; 

['Relative to fundamental < ']; 

[ 'Harmonic Number' ] ; 

[ 'Normalized frequency' ] ; 

I t • 

t 

['Srd harmonic power sampled and quantized = '] 
[textlO num2str(LPP(2) ) • dBr']; 

['Signal freq = ' num2str(fo) ', clock freq = 
[textll num2str (fc) ] ; 

['Srd harmonic power quantized only = ']; 

[textlS num2str (pp(2) ) ' dBr']; 



] ; 



lOS 



f2 = zeros (Maxhar, N) ; 
disp ('calculating the q') 
for m = 1:5, 

A3 = m*tc*wc/2; 
sinc(m) = sin(A3)/A3; 
for n = 1:100, 

freq = m*fc - (2*(n-l)+l)* fo; 
if abs(freq/fo) <= Maxhar 

harno = fix(abs(freq/fo)/2)+l; 
for k = 1:N, 

since (harno) = sinc(m) ; 
hhno (harno) = freq/fo; 

s(n) = cos( (2*(n-l)+l) *asin( (k-l)/(N) ) ) ; 
hno(n) = 2*(n-l)+l; 

q(n,k) = Al* A2/ (2* (n-1) +1) * sine (m) *s (n) ; 
f = f 2 (harno, k) .* f 2 (harno, k); 
f2(harno,k) = (f + q(n,k) .*q(n,k) ) .^0.5; 

end, 

end, 

end, 

end, 

disp ( 'calculating QK') 
qmnk = zeros (l:Maxhar) ; 
for k = 1:N, 

qmnk = qmnk + f2(:,k)'; 
end, 

PP = (pnk'+qmnk)/ (pnk(l) +qmnk(l)); 

LPP = 10*logl0(PP .* PP) ; 
normfreq = hno*fo; 

plot(normfreq(l:harp/2+l) , LPP( 1: harp/2+1) , 'o') 

title ( [num2str (bits) ' bit amplitude quantized sampled 

signal ' ] ) 

xlabel (textSl) 

ylabel (text?) 

grid 

text (0.2, 0.75, textlO, 'sc') 
text (0.2, 0.8, textll, 'sc') 
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text (0.2, 0.7, textl3, 'sc') 
meta drl 



% FOR PRINTING THE PLOT 
%!gpp drl 



%FOR CONVERTING THE PLOT COMPATIBLE WITH WORDPERFECT 
!gpp drl /dhpgl /f drl. pic 
disp('pic file is DRl.pic') 



filename = 'SAM.dat'; 

disp ( 'Writing data to file SAM.dat') 

freqt = ['Normalized input signal frequency = ' ]; 

freqt = [freqt num2str(fo) '\n\n\n']; 

pw = ['Clock pulsewidth to period ratio = ' num2str(tc) 

noofbits = [ 'Number of bits of quantization = ' ] ; 

noofbits = [noofbits num2str (bits) '\n']; 

fprintf (filename, noofbits); 

fprintf ( filename , pw) ; 

fprintf ( filename, freqt) ; 



'\n'] 



ptextl = 
ptextl = 
fprintf ( 
ptext2 = 
fprintf ( 
fprintf ( 
ptext3 = 
fprintf ( 
ptext4 = 
fprintf ( 



[ptextl 

filename, 

'Harmonic 

filename, 

filename, 

I 

filename , 

t 



Power level relative to' 

' fundamental\n' ) ; 
ptextl) ; 

Sampled Quantized 

ptext2) ; 

' number 

ptext3 ) 



quantized 



\n' 



only\n' ) 



signal 



signal \n\' 



dBr 



dBr 



\n\n\n' 



filename, ptext4) 



% MAXIMUM OF 19 HARMONICS FOR FILE DATA 
if harp/2 >= 10, 
kmax = 10; 

else 

kmax = harp/ 2; 
end, 

ptextS = '%4.0f %18.2f %16.2f \n' 

for k = irkmax 

fprintf (filename, ptextS, hno(k) ,LPP(k) ,pp(k) ) 
end. 
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